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A MATKEMATICAL TBEORT OF PLASTICITI BASED 


ON THE CONCEPT OF SLIP 
By S, B, Batdorf and Bernard Budiansky 


SUMMART 


A theory of plasticity ba.sed on the concept of slip is proposed 
for the relationship between stress and strain for initially isotropic 
materials in the strain-hardening range. As treated in the present 
paper, the theory is an extension to polyaxial stress conditions of 
the conventional uniaxial stress— strain relation, and time-dependent 
effects, such as creep and stress relaxation, are not considered. 

Previously existing theories can in general be classified into 
two groups, often called, theories of plastic deformation and theories 
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that for continual loading the state of strain is imlquely determined 
by the state of stress. Flow theories are based on the assumption 
that the increment of strain is unlq^uely determined by the existing 
stress and the increment of stress. Theories of both types make 
additional assumptions concerning the orientation of the principal 
axes of plastic strain or of Increment of plastic strain and concerning 
the nature of loading and unloading. 

The theory proposed is of neither the flow nor the deformation 
type . Its formulation vas guided mainly by physical considerations 
with respect to the assumed mechanism of plastic deformation, namely 
slip. The new theory makes no use of the assunqjtions just mentioned 
in connection with flow and deformation theories; in fact, according 
to this theory, all these asaun^tlona are wrong. 

The new theory is based on the assumption that slip in any direc- 
tion aJLong parallel planes of any particular orientation in the material 
gives rise to a plastic shear strain which depends only on the history 
of the corresponding component of shear stress. The relation between 
this plastic shear strain and the corresponding stress can be determined 
from the tensile stress— strain curve. The plastic strain due to any 
system of applied stresses is then found by considering the history of 
the con5>anent in each direction of the shear stress on each plane of 
the material, finding the corresponding plastic shear strain, trans- 
forming this plastic shear strain into plastic strains in some fixed 
system of coordinates, and summing over all slip directions and slip 
plane orientations. A semigraphical method is given for computing in 
this manner the strains associated with a given stress state. 
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The theory is shown to give results in tetter agreement than 
previously existing theories with data obtained in an experiment in 
which a cylinder was compressed into the plastic range and then twisted 
at constant compressive strain. The theory therefore appears promising, 
but more theoretical and experimental investigations of the character- 
istic differences between the new theory and previously existing 
theories are required before the evidence for or against the new theory 
can be considered conclusive. 


INTROrfCJCTION 


A new theory is proposed for the relationship between stress and 
strain for initially isotropic metals in the strain-liardening range. 

The theory seeks to predict the plastic strain that would resiilt from 
the application in any sequence of arbitrary combinations of stress 
and requires for this purpose only a knowledge of the uniaxial stress- 
strain relation for the material. 

This theory is essentially quite different from the many theories 
of plasticity that have previously been proposed. It has been recognized 
(reference l) that existing theories of plasticity may in general be 
classified into two types: deformation theories and flow theories. 

Both types of theory assume that the relation between stress and strain 
is governed by one law during loading and by another law (the elastic 
law in incremental form) diuring unloading. Loading and unloading are 
defined in terms of the increase or decrease of some fotationally 
Invariant function of the stresses, those most frequently used being 
the maxiinum and the octaiiedral shear stresses. The two types of theory 
differ chiefly in that deformation theories propose laws giving plastic 
strain in terms of applied stress; whereas flow theories propose laws 
giving increments of plastic strain in terms of applied stress and 
Increments of applied stress. Thus, the previous loading history is 
assumed to have no effect on strain in deformation theories or on 
Increment of strain in flow theories. These basic assumptions appear 
to be in the nature of plausible postulates made to facilitate the 
mathematical formulation and subsequent analysis of the various 
theories rather than ph.ysical.ly evident or experimental 1 y verified 
facta. Moreover, it is characteristic of both flow and deformation 
theories that the laws are proposed by means of equations having forma 
that are more or less arbitrary, which therefore introduce additional 
assumptions of uncertain validity. For example, in theories of plastic 
flow, the equations proposed are such that the principal axes of stress 
and of increment of plastic strain are assumed to coincide. On the 
other hand, deformation theories assume, coincidence of the principal 
axes of stress *and plastic strain. Consequently, although canslderable 
difference of opinion exists concerning the relative merits of flow and 
deformation theories, it is entirely possible that neither type of 
theory is correct. 
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The present theory is of neither the flow nor the deformation type. 
INilike most theories proposed to date, its formulation was guided mainly 
hy physical considerations with respect to the mechanism of plastic 
deformation. The present theoiy involves none of the assumptions Just 
discussed in cannection with flow and deformation theories; Indeed, 
according to this theory, all these assunqjtlone are wrong. 

This paper presents the essential features of the new theory and 
seeks to establish its claim to serious consideration hy showing its 
marked superiority to previously existing theories in accoxmting for 
the data obtained in one of the more crucial types of plasticity 
experiment . 


SYMBOLS 


solid angle, used to describe orientation of slip 
planes, steradlans 

angular coordinate giving direction of slip, radians 
direct stress 

direct stress at Tdiich plastic defonaatlon begins 
shear stress 

shear stress at which plastic deformation begins 
total direct strain 
plastic strain 
7 total shear strain 

7 " plastic shear strain 

X, y, z rectangular coordinates; used also as subscripts in 

cannection with stress and strain to denote 
particular components of these quantities 

1, 2 direction of normal to slip plane and direction of • 

slip, respectively 

characteristic shear function for material, giving 
plastic shear strain per steradian of slip-plane 
orientation per radian of slip direction as a 
fmction of shear stress 

1 x 1 j ^yl^*«»^z2 cosines of angles between x and 1, y and 1, ... 

z and 2 directions 


F, f(tj^) 

4 



I 


mm 


liliilippp 
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coefficient of nth. term In series expansion for F 


/g \ 


function giving variation of plastic strain with 

applied stress in uniaxial compression or tension, 
corresponding to nth term of series expansion for P 


n^, Dg 


mit vectors In 1 and 2 direct lens, respect lvel 7 


1, J, k 


unit vectors in the i, y, and 2 dlrectlans, 
respectively 


Q. 


unit vectors having dlrectiona indicated in figure 8 


a, 0 


polar coordinates used to specify orientation of slip 
plane (see fig. 8 ) 


ASSOMPTIONS AND UKCTATiaNS 


The uniaxial stresEKStraln hehavlor assumed in the present analysis 
is the subtly idealized one usually adopted for engineering purposes. 
The stress— strain curve for first loading is represented hy OABC in 
figure 1. The straigjit part of this curve QA is called the elastic 
range, and the point A is called the elastic limit . If the load is 
removed "before the stress exceeds the elastic limit, the material returns 
to its original state at 0, and there is no permanent defoimation. If, 
on the other hand, the stress is increased to the value at point B and 
then is removed, the stress-strain relationship during the unloading 
process is that given by the dashed line BO* . The line BO* has the 
same slope as the line QA ,so that during the unloading the material 
behaves elastically except that it has undergone permanent plastic strain 
given by 00*. Upon subsequent loading, the point representing the 
state of the material stays on the dashed line until it reaches B and 
then, as the load further increases, continues along the original path 
toward C. Additional unloadings and loadings follow the same pattern 
of unloading along a line parallel to OA and loading along the same 
line until the stress-strain curve of the original material is reached. 

At any stage of the loading history Just described, the strain may be 
.considered to be the sum of an elastic part related to the stress by 
Hooke *s law and a plastic part determined by the hipest value of the 
stress occurring up to the time under consideration. 


The law of plastic deformation to be postulated is in accordance 
with the uniaxial behavior Just described. In its original unstressed 
state the material is assumed to have identical stress-strain curves in 
tension and compression. The effect of reversal of loading upon the 
plastic behavior la not treated in this paper other than to state in 
passing that the theory is flexible enou^ to be adapted to materials 
with a Bauschinger effect. As in the uniaxial case, time— dependent 
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effects, such as creep, elastic recovery, and the effect on the stress- 
strain curve of rate of loading, are ne^ected. 


THTT ORV 


Physical ccnsideratione .— Slip lines are cammoniy ohserved in metals 
undergoing plastic deformation. On the assumption that, at least for 
some metals, slip is the principal mechanism of plastic deformation, the 
present theory was formulated in such a manner as to incorporate the 
main features of this phenomenon. Those features are now reviewed. 

Vfhen a metal crystal is subjected to a rnnall shear stress, it undei^ 
goes elastic shear deformation. (See fig. 2.) As the stress increases 
beyond a certain limiting value, small blocks of the crystal commence to 
slide with respect to each other el nng crystallographic planes called 
slip planes. As the displacements increase, the shear stress required 
to produce further deformation also increases. The total shear displace- 
ment at any time is the sum of the elastic shear in the crystal blocks 
and the shear displacement due to slip between nel^borlng blocks. 

When the stress is removed, the elastic componant vanishes, but the 
displacements due to slip remain as a plastic strain. The relationship 
between applied shear stress and total shear strain for repeated loading 
and unloading is similar to that described for tmiaxial tension end 
illustrated in figure 1. The shear stress required to produce slip is 
found to be substantially Independent of the normal stress and slip 
occurs only along certain preferred directions in the slip plane. (See 
reference 2.) 

Accordingly, in the theory, the following assunqjtlonB are made: Any 

strain is the sum of an elastic strain and a plastic strain; plastic 
strain is conq)oeed solely of shear deformations due to slip; these shear 
deformations are uninfluenced by normal pressures; the plastic shear 
deformation resulting from slip in any direction in a plane of any 
given orientation depends only upon the history of the component in the 
direction of slip of the shear stress on this plane; and the total 
plastic strain is sinqily the sum of all the slip deformations that have 
occirrred. If the stress history is such that reverse slip does not 
occur, the dependence on the aforementioned previous history becomes 
eiiiq)ly a dependence on the highest previous value of the sheai>-stres3 
canq>onent in question. 

Mathematical formulation of the theory .— Althou^ a macroscopically 
isotropic metal is actually an aggregate of tiny crystals, it is 
customary and convenient to foimulate its •stress-strain relatione as 
though it were a continuum. Accordingly, instead of considering a sman, 
slip along each of a large number of discrete planes, the theory 
contemplates an infinitesimal plastic shear strain associated with 
each infinitesimal fraction of the contlnuiom comprising all possible 
planes . 
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It is convenient to represent the orientation of a particular plane 
hy the coordinates of the point at vhich it would he tangent to a hemi- 
sphere. The radius to the point of tangency is the normal to the plane. 
The axis 1 shown in figure 3 ie tahen in the direction of the normal to 
the plane, and the axis 2 denotes a particular direction of slip in the 
plane. An infinitesimal hand of planes may he represented hy the normals 
included in the solid angle dC, and an infinitesimal hand of slip 
directions may he represented by d3. The theory postulates that the 
slip along the planes dfl, in the increment dp of slip directions, 
produces an infinitesimal plastic shear strain d/ip" associated with 
the 1,2 axes that is given hy 




where F is a function depending only on the history of T^p, the 
shear stress in the 2 direction on the plane perpendicular xo the 
1-axls. 


The contribution of this infinitesimal shear strain to the strains 
in the standard i— , y— , and z-axes is readily written in terms of the 
direction cosines of the 1-axis and 2-axls. Thus, 


^1" * ^xl^i 2 ^^12” 


d€y" = *^’12" 


^zl^z 2 ^^12" 


d7 " 
3cy 


(^xl^yS Vl^x2)^^12" 




= (^il^z 2 ^zl^i^ ^^12” 


^^yz" 




In order to find the total plastic strains in the standard axes, the 
effects of the plastic sheai>-etrain Increments must he integrated over all 
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directions and all planes. Thus, 


=l" ■ / / , “ W 

H -^/2 



(3) 


•where H denotes the entire hemisphere. {The slip-direction increment 
is integrated through l 80 ° to include all possible directions only once. ) 

If the applied stresses are giTen in the standard coordinate system 
as 0 ^, Cy, . . . Ty 2 * shear stress T -,p is given by 


■^12 " + hi^z2°z + (^ilVs Vl^xa^ry 


0xl^z2 ^z1^x2)tiz + (Vl^z2 ^zl^yaVyz 


In order to e'valuate ■the total plastic strains as given by 
equations (3), analytical expressions for the direction cosines must 
be used. Appendix A contains such expressions in a system of spherical 
coordinates. 

Characteristic shear function .— It follows from the previous 
discussion concerning slip that, if T-,p is gradually increased, the 

function F must remain zero until a limiting value Tj^ is reached. 
This limiting value is evidently equal to one-half the elastic limit Oj^ 
in pure tension (or compression) since in a tension test the maximum 
shear stress is one-half the applied uniaxial stress. Beyond this 
value, F varies with ^ ^ ™ann©r characteristic of the material 


< 
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and is referred to hereinafter as the "characteristic shear function." 
and the corresponding curve, as the "characteristic shear curve." 

(See fig. 4. ) 

The characteristic shear curve can he determined from the stress- 
strain curve for the material in the following manner. The characteristic 
shear function F is expanded into the series 




•which is then substituted into equation ( 3 ). The resulting plastic 
tensile strain can he shown to he given h 7 



where the functions do not depend on the material under considera- 

tion. A set of these functions for the case N = 5 has been evaluated 
hj the method described in appendix B over a range of stress extending 
up to 1.8 times the elastic limit stress. The results are tabulated in 
table 1 and are shown graphically in figure 5» 

The determination of the characteristic shear function corresponding 
to a particular stress— strain curve by means of the g-functions is 
quite slnple. Equation (6) is written for 0 / 0 ^ ®q.ual to 1.1, 1.25, 

1.4, 1.6, and 1.8; table 1 is used to obtain the numerical values of 
the g-functions at the corresponding strain in each case. This procedure 
gives five linear algebraic equations for determining the five unknown 
coefficients a^^, ag, . . . a^. These values of the coefficients are 
then substituted into equation ( 5 ). If the stress-strain curve is not 
known to a stress as hi^ as 1.8 times the elastic limit stress, either 
fewer equations and fewer g-functicne or interpolated values of the 
g-functions must bo used. 

Application of theory .— Once the characteristic shear curve has 
been obtained, it is only necessary to apply equations ( 3 ) and (4) to 
obtain the plastic strains resulting from a given sequence of stresses, 
^e tj 2 used to evaluate F is, as indicated previously, the largest 
value of this stress which has occurred for each particular choice of 
directions 1 and 2 during the history of the loading. Because serious 
difficulty is Involved in evaluating the integrals analytically, the 
use of approximate methods is advisable. Appendix C describes a semi- 
graphical method of evaluating the required Integrals for any plane— 
stress condition. The method provides a stral^tforward procedure for 
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the calculation of strains with accuracy sufficient for engineering 
purposes. 

COl-IPAEISON WI'Ei EXPSRIMEaSTAL DATA AND DISCUSSION 


In mc5st of the experiments that have heen performed to test theories 
of plasticity, the ratios and directions of the principal stresses hare 
heen constant throu^out the loading. In general, such experiments fail 
to distinguish sharply between the various theories, and in particular do 
not provide any basis for choosing between flow and deformation theories 
(references 3 and 4). 

A more crucial type of experiment was perfoirmed by Eoger W. Peters 
and Norris F. Dow at the Langley Laboratory of the NACA. In this experi- 
ment, a thin aluminum-alloy cylinder was compressed into the plastic 
range and, with the conpressive strain held constant, was then twisted. 
During the twisting process, the compressive stress decreased. After a 
shear stress of nearly 12 ksi had been applied, the shear stress and 
then the compressive stress were removed. 

In order to facilitate con^jarison between the experimental data 
and the predictions of various plasticity theories, the experimental 
stresses were adopted as the prescribed conditions and the corresponding 
strains were computed. The stress combinations applied during twisting 
and the resulting plastic compressive strain are shown as solid curves 
in figure 6. It is evident from the figure that plastic action continued 
throughout the twisting process. This behavior is contrary to what 
would be espected on the basis of conventional theories. Consider, for 
example, the dashed curve representing constant octahedral shear stress 
In figure 6. The region below this curve corresponds to octahedral shear 
stresses lower than that at the beginning of twisting, and the region 
above the curve corresponds to hi^er values of the octahedral shear 
stress. Comparison of the stress-history curve and the curve of constant 
octahedral shear stress shows that the octahedral shear stress was 
decreasing during the early stages of twisting and did not Increase 
above the hipest previoua value until the applied shear stress exceeded 
8.6 ksi. On the assumption that after unloading, plastic action is not 
resumed until the hipest previous value of the loading function is 
exceeded, octahedral— shear theories would give zero plastic strain up 
to a shear stress in excess of 8.6 ksi, and similarly maximum-shear— 
stress theories would predict zero plastic strain up to a shear stress 
of about 5 ksi. The present theory, on the other hand, predicts plastic 
action throu^out the twisting process, because even when the maximum 
shear was decreasing in magnitude, its orientation was being changed and 
new and unhardened planes were being subjected to high shear stresses. 

Figure 7 shows that the present theory is not only qualitatively 
but also quantitatively in rather satisfactory agreement with the test 
data, in any event in much better agreement than conventional theories. 
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The toteO. plastic ccm^jresBl've strains, ccai 5 >uted "by the present theory 
at shear stresses of 6,3 ksi and 11.8 ksi, were helow the measured 
values by 6 percent and 9 percent, respectively. As the figure indicates, 
the present theory predicts an almst constant total con^jressive strain 
in agreement with experiment. The octabedral-ehear theories would 
predict the elastic relation up to an applied shear stress of 8.6 ksi 
a-ni^ the maximuift-ehear theories, to an applied shear stress of about 5 ksi. 
The predictions of these theories at hi^er applied shear stresses have 
not been plotted because there appears to be no general agreement as to 
the proper way to apply them in the case of second loading. 

The possibility exists, of course, that new flow or deformation 
theories might be devised using a new stress invariant that would increase 
continually during the loading sequence accompanying twisting in the 
experiment Just described, and it mi^t at first be suspected that the 
theory of the present paper is perhaps a flow or deformation theory 
en^loying such a loading criterion. The fact that this is not the case 
can be seen by considering an experiment in which a state of pure tension 
in the plastic range is maintained approximately constant in ma g n itude 
■vdiile being rotated in direction, (See reference 4.) According to both 
flow and deformation theories, if the magnitude of the pure tension 
decreases slightly during rotation, no plastic deformation can occur. 
According to the present theory, however, new and unhardened planes 
acquire large- shear stresses during the rotation, and plastic deformation 
occiars. Moreover, by definition (reference 5 ), deformation and flow 
theories predict the Increments of the strains for continual loading on 
the basis of only the corresponding increments of stress and the state 
of stress immediately prior to application of the stress increments; 
whereas the present theory requires for such a computation a knowledge of 
the entire previous stress history so that the extent and nature of the 
hardening of each plane can be deteimlned. 


CONCLDDING REMARKS 


Previously existing mathematical theories of plasticity can generally 
be classified into two types, often called flow and deformation theories. 
Both types of theory are constructed largely on the basis of mathematical 
considerations mifl involve a number of arbitrary assumptions of uncertain 
validity. The theory proposed herein is of neither the flow nor the 
deformation type but constitutes an entirely different approach to the 
problem of plastic behavior, its formulation being guided mainly by 
physical considerations with respect to the assumed mechanism of 
plastic deformation, namely slip. 

Test data obtained in one of the more crucial types of plasticity 
experiment are found to be in decidedly better agreement with the present 
theory than with previously existing theories. The experimental data 
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hovever stre not very , extensive . Before the new theory can be regarded 
as established, additional e3i3)eriiiiental data of the same general type 
must be obtained, the other characteristic differences between the new 
theory and other theories of plasticity must be explored both theoreti- 
cally and experimentally, and the range of stresses and strains within 
f^lch the theory is reasonably accurate must be determined. 


Langley Aeronautical Laboratoiy 

National Advisory Committee for Aeronautics 

Langley Air Force Base, Ya., February l4, 19^9 
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APPENDIX A 


DERIVATION OF TRANSFORMATION EQUATIONS IN SPHERICAL COORDINATES 


In order to find the shear stress In a given slip plane due to a 
set of applied stresses as well as to perform the Integrations required 
to determine the corresponding plastic strain, analytical expressions 
for certain direction cosines used In the analysis should he obtained. 
Such ea^resslons are now obtained for a polar-coordinate system In 
which the z-ails Is chosen as the polar axle. Because the plasticity 
theory under consideration makes use of only two directions — namely, 
the direction of slip within the slip plane and the normal to the 
plane — only the cosines of these directions are derived. The corre- 
sponding unit vectors n^ and n^, as well as the unit vectors T, 

J, k, m, and q are identified in figure 8. 


The direction cosines 1. 


of the unit vectors 
Thus, 


■xl> 
n^ and 112 


Vl* * ' 

In the 


1^2 slDjply the components 

y, and z directions. 


and 




112 = 1^x2 + J^y2 ^'^z2 


(Al) 


This resolution of the vectors may be brou^t about in a succession of 
easy stages as follows: 


From figure 8 


n-j_ = k sin 0 ->■ p cos 0 
p = 1 sin a + J cos a 


(A2) 

(A3) 


Prom equations (A2) and (A3) 

n^^ = T sin a cos 0 + J cos a cos 0 + k sin 0 
Similarly, from figure 8 

IS 2 = cos 3 — ^ sin p 


(a4) 


(A5) 
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m = k cos 0 — p sin 0 

(A6) 

q = —1 cos a + J sin a 

(A7) 


From equations (A3), (A5), (a 6), and (A?) 

n^ = l(cos a sin ^ — sin a cos 3 sin 0) 

+ j(— Bin a sin 3 — cos a cos 3 Bln 0) 

+ k cos 3 cos 0 (a8) 

From equations (A1), (Alt), and (a 8) the eipressions for the direction 
cosines are identified as follows: 


= sin a cos 0 


2yl = cos a cos 0 

> 

^zl = sin 0 




- COB a sin 3 — sin a cos 3 sin 0 I 

Zy 2 = —sin a sin 3 — cos a cos 3 sin 0 ? 

1^2 - P ^ 



(A9) 


(AlO) 


Suhstitutlon of the values for 2^.^^ * 1 » 1^2 equations (A9) 

and (AIO) in equations (4) and (2) of the text leads, after some manipula- 
tion, to the following resiilts; The component in the direction ng of 
the shear stress in the plane normal to n]^ due to the application of 
stresses o^, Cy, ... Tyj. is given hy 


T .|p = i a^(sln 2a sin 3 cos 0 — sln^a cos 3 siJi 20) 

+ i a^(-ein 2a sin 3 cos 0 — cos^a cos 3 sin 20) 

+ 20 COS 3) 

+ T^(cos 2a sin 3 cos 0 — ^ sin 2a cos 3 sin 2J^ 

+ T2 ^(cos a sin 3 sin 0 + sin a cos 3 cos 20) 

+ T (—sin a sin 3 slJi 0 + cos a cos 3 cos 20) (All) 

yz 
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The increments of strain resulting from the increments of shear deformation 
in the slip plane d 7 i 2 " given hy 




= P dy.jp"(ain 2 a sin 3 cos 0 — sin^a cos P sin 20 ) 


de/ 


= ^ d 7Tp "(— Bin 2 a sin p cos 0 — cos^a cos P sin 20 ) 
= i d72_2"(ain 20 cos p) 

= d722"C°°® ^ 0 ~ 2 ^ 




" 
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(A12) 


= d7i2"(cos a sin p sin 0 + sin a cos p cos 20) 


d7y2" = d722"(“fl^ ® sin p sin 0 + cos a cos p cos 20) 
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APPMDIZ B 


DElEERMIMTiail OF GHAEACTERISTIC SHEAH FUNCTICBf 


The characteristic shear function for a material may te conveniently 
determined from the uniaxial stress— strain curve for the material. If 
the 2-€txls is the ails al.ong vhich stress is applied (see fig, 8), 
equations (3), (A9), and (AlO) ccmblne to give for the corresponding 
plastic strain 


pjt/2 pjr/2 

= / d.0 d3 cos 3 sin 20(2n cos 0) (Bl) 

'^0 '^- n /2 2 


This integration is perfonned over the hemisphere in front of the x,j plane 
and, since axial symmetry exists, dfl is taken equal to 2it cos 0 d0. 


The characteristic shear function F is expanded into the series 

K 


^’(^12) = ®n 

22^1 


&-*)■ 


(B2) 


for 


'12 ^ "’’L* 

Consequently, idien the value of F from equation (B2) is suhstltuted 
into equation (Bl) the integration should he carried out only over the 
regions "where T-|p ^ or 


For T-|p < Tj^, the characteristic shear function is zero. 


a sin 20 cos 3 > o_ 


(B3) 


because for uniaxial tension 


'12 


= i a sin 20 cos 3 


(BU) 


and 




= 2 


(B5) 
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From the use of equations (B 2 ) to equation (Bl) hecomes 


p jt— t p ^ N ~j 

~ ^ ^ ^ P “ l)^sin t cos cos ^ dp dt 

->t ^ -p itel I 


•where 


t « 


sin ^ — 


1 


8 sin t 



The integration -with respect to g in equation (b 6 ) can he caixied out = 
analytically and leads to -the following result: 


s - 5 


where 


oK-sin”^^^ 

e^(s) = :t / 

J sin”^( i) 

\S/ 
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and 




= cos 


^ sin t|^ s sin t(p + sin p cos *^) — sin 


= cos sin tjj s^sin^t sin p(cos^ + 2) 

— 8 sin t(p + sin p cos p) + sin pj 

= cos sin tjs sin t(p + sin p cos s^sin^t + 0 

— s^sln^t sin p(2 + cos^) — ^ sin ^ 


y (B9) 


q^ = cos sin tjs^sin^t sin p(2 + cos^)^^ s^ain^t + ^ 
— s sin t(*^ + sin p cos s^sin^t + 2^ + j sin 

q^ = cos sin t|s sin tC? + sin cos s^sinH 


15 ? 




s^sin^t + ^ - s^sin^t sin p(2 + cos^)^^ s^sin^t + 


+ sin Plfr s^sin^t + 


J 


The integrations in equations (b 8) have heen carried out numerically 
for the case N = 5 for values of equal to 1.1, 1.25, 1.4, 

1.6, and 1.8. The integrations were performed hy finding the area under 
fifth-degree curves passing through the first six and last six points 
corresponding to 11 values of 0- at each stress level, seven decimal 
places being retained during the computations. The results obtained are 
given in table 1 and are shovn graphically in figure 5* 




NACATN No. I 87 I 


The e(iuatlcai 




(BIO) 


should glre a very close approitmatlom. to the characteristic shear 
function for the material when the values a^^, a^, . . . a^ are chosen 
to satisfy equation (B7) for the case N = 5 five stress levels 

for •vdiich the niunerical Integrations were performed. Except for the 
errors In the integrations, the characteristic shear function so found 
would lead to a stress— strain curve fitting the one from which it was 
derived over the entire elastic range and at the five points in the 
plastic range already Indicated. If desired, of course, less accurate 
characteristic shear functions can he found with less confutation hy 
using fewer g-functions and hy fitting the stress-strain curve at 
fewer points. 
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APPENDIX C 


A METHOD OF DETERMINING PLASTIC STRAINS ASSOCIATED 
WZCH PLANE STRESS 

As indicated in equations (3) of the text, the theory gives plastic 
strains in terms of triple integrations. Because these integrations would 
prove quite formidable if treated analytically, consideration has been 
given to approximate methods. The method now to be described is a semi — 
graphical one designed to facilitate the con^iutaticm of plastic strains 
associated with plane stress. 

tvia system of plane stress is assumed to act in the x,y plane, 
ymft the hemisphere on idiich the directions of the slip— plane normals 
are measured 1s taken above the i,z plane. This hemisphere is considered 
to be divided by lines of constant a and constant 0 into elements of 
equal area, the projections of which on the x,y plane are illustrated in 
figure 9. The elements of area on the front half of the hemisphere are 
denoted by a system of subscripts as indicated in the figure. The mean 
coordinates each element of area were calculated and are 

listed in table 2. 

Tn order to find the shear stress 7^2 ^ plane tangent to the 

hemisphere at the point equation (All) is written in the form 


( 7 ^ = fK-Vi 3 - Bnn COB p) 

\ L/mn I* 


- -^Ajjm sin 3 + Cn^i cos 3) 


+ Bln 3 - Enm cos 3 ) 


where 


t 


Anm = Bin 20ni cos 0^ 
Bujp = sln^Om sin 20^ 


(C3) 
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Cnn = cos^c^j^ sin 20^ 

(C4) 



= cos 2oni cos 0^ 

(C5) 



= I sin 2oni sin 20^ 

(C6) 



The numerical values of these coefficients are given in table 3. 

The contributions of the plastic shear deformation arising from the 
slip in planes tangent to the part of the hemisphere in the m,nth element 
of area to the plastic strains measured in the system of coordinates 
used to specify the applied stresses are given by a similar reformulation 
of equations (A12): 


C^x ^mnOViia P "" ®ran P) 

(C7) 


(d€y”}mn = “ ^*^^12")mnOWi P + ^mn cos 3) 

(C8) 

■ 

(*^^3cy")mn ~ (d7ip")Tim6^mn P ~ ®nm 

(C9) 

/ 


Equations (C7) to (C9) make use of the assunq)tion, later to be employed 
also in applying equation (Cl), that all planes tangent to the part of 
the hemisphere in any rectangular element of area may for computational 
purposes be regarded as concentrated at the mean coordinates of the 
element. 

The values of the shear- stress component a-'t point (m,n) in the 

directions p = 0° and 3 = 90° are computed by using equations (Cl). 
These components are then plotted at ri^t angles from a common origin, 
and a circle is drawn throu^ the origin and two ends. (See fig. 10 (a).) 
The diameter drawn from the origin then represents the magnitude and 
direction of the shear stress in the m,nth plane, and the chords drawn 
from the origin in othdr directions represent the conq>onents of shear 
stress in those directions. 

A circle is then drawn about the origin with a radius corresponding 
to the elastic— limit shear stress, as in figure 10 (b). According to the 
theory, shear deformations are associated with each direction in which 
the shear-stress component exceeds the elastic— limit shear stress, the 
amount of shear deformation in any direction depending on the amount by 
which the corresponding shear— stress component exceeds the elastic— limit 
shear stress. 


In summing the contributions in the different directions, only 
discrete directions separated by 9°» s-s indicated by radial lines in 
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figure lOCt), are considered, 
directions Is given by 


Hie shear strain In each of these 





and Is found by using the sheaj>-etrees component In that direction and 
the characteristic shear curve for the material. (See fig. 4.) Eefer— 
ence to equations (C7) to (C9) shows that the products of this strain 
with sin p and with cos P must then be found. The summation over 
the various directions of slip (direction 2) may then be carried out; 
this procedure coirresponds to performing the integration over P in 
equatione (3). 

The operations of finding the plastic shear deformation in each 
direction and multiplying the result by sin P or by cos p may 
conveniently be carried out in a single step by the use of a master 
chart of the type shown in figure 11, Along each radial direction 
considered are plotted values of F cos P for the material being used. 
This chart Is simply fitted over figure 10(a), -vdiich represents the 
stress state, and the values of F cos P are read off directly (see 
fig. 12) and added in order* to perform the integration over p. The 
value of this sum is then multiplied by the appropriate coefficient, 
determined from equation (C7), (c8), or (C9), and table 3, The 
values of F sin p may be obtained in the same manner by rotating 
the master chart throu^ 90° • A different work sheet is used for 
each point (m,n) at which the shear stress on the tangent plane 
causes plastic shear deformation. For plane— stress systems symmetry 
exists between the front and rear halves of the hemisphere; thus, 
the total plastic strains ^y”, and are found by adding 

up the corresponding strains computed on each sheet associated with 
the front half of the hemisphere and by multiplying by 2. 

At first si^t this method of integration appears quite formidable 
because the front of the hemisphere is divided into 200 elements, which 
would appear to require 200 work sheets. In pi^ctice, however, it is 
found that ordinarily many of the elements (m,n) make no contribution 
to the plastic deformation. In addition, the symmetry properties of the 
stress state considered are often such that only a quarter or sometimes 
only an eighth of the hemisphere needs to be considered. 
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TABLE 1 

CALCULATED g-^ALUES 


«?Al 

81 

82 

«3 

84 

«5 

1.10 

1.25 

1.14-0 

1.60 

1.80 

0.03111^86 

.1671619 

.37^^986 

.72III448 

1.1965026 

0.0020860 

.0281396 

.iooia47 

.29I48256 

.6148685 

0.0001560 

.0053115 

.0308275 

.1351^768 

.3883228 

0.0000137 

.001071a 

.0099147 

.0651761 

.2I487629 


0.0002206 

.0033785 

.0336614 

.1656213 





NACA TN No. 187 I 


TABLE 2 

MEAN COGRDIKATES OF ELEMEBTS OF AHEA 

[see fig. 9; = -OjJ 


m^n 

(deg) 

(4eg) 

1 

1^.5 

2.9 

2 

13.5 

8.6 

3 

22.5 

14.5 

k 

31.5 

20.5 

5 

40.5 

26.8 

6 

49.5 

33.4 

7 

58.5 

40.6 

8 

67.5 

1 * 8.8 

9 

76.5 

58.6 

10 

85.5 

“ 72.8 


^Taken one— third of •wa 7 from houndary 


of triangular element to pole. 


























































Figure 2- Slip in single crystal (schematic) 
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Figure 5.- Variations of tensile plastic strain with 
applied tensile stress corresponding to various 
terms in series expansion for characteristic 
shear curve. 
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Rgure 6.~ Compressive stress and plastic compressive 
strain measured during twist of cylinder at 
constant total compressive strain. 




Figure 7- Compressive stress-strain relations for cylinder subjected first 
to compression (O/i), and then to torsion at constant compressive 
strain (AB), followed by removal of torsion and then compression. 
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Figure 8.- Coordinate system used for analytical 
representation of direction cosines. 




Figure 9.- System for classifying slip planes by means 
of plane orientation (schematic). 
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(a) Components of sheor stress in various 
directions of slip. 



(b) Slip directions ore those in crescent. 

Figure 10.' Graphical construction used in determining sheor 
stress components cousing slip. 
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F(ij,)co»;9 « 10* 

Figure II.- Master chart for rapid evaluation of F(r,^)cos/9 and F(r|^)sin^. 
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Figur* 12. ~ Matt«r chart fitted over cheor-etreee diogram 
(fig. 10(a)) to permit direct reoding of F(r, 2 )coe /9 
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